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ABSTRACT 
The Baryonic Tully-Fisher Relation (BTFR) is an empirical relation between baryonic mass and 
rotation velocity in disk galaxies. It provides tests of galaxy formation models in ACDM and of 
alternative theories like MOND. Observations of gas rich galaxies provide a measure of the slope 
and normalization of the BTFR that is more accurate (if less precise) than that provided by star 
dominated spirals, as their masses are insensitive to the details of stellar population modeling. Recent 
independent data for such galaxies are consistent with Mb — AVf with A = 47± 6 M© km~ s"*. This 

is equivalent to MOND with oq = 1.3±0.3 As^^. The scatter in the data is consistent with being due 
entirely to observational uncertainties. It is unclear why the physics of galaxy formation in ACDM 
happens to pick out the relation predicted by MOND. We introduce a feedback efficacy parameter £ 
to relate halo properties to those of the galaxies they host. S correlates with star formation rate and 
gas fraction in the sense that galaxies that have experienced the least star formation have been most 
impacted by feedback. 

Subject headings: galaxies: kinematics and dynamics — galaxies: dwarf — galaxies: irregular — 
galaxies: spiral 



1. INTRODUCTION 



The Tully-Fisher relation (jTuUv fc Fished Il977| ) is 
one of the strongest empirical correlations in extra- 
galactic astronomy. I n addition to b eing a useful 
distance indicator fe.g.. iSakai et al.l 120000 . the physical 
basis of this relation is of fundamental importance. 
In the context of ACDM, the Tully-Fisher relation 
places tight constraints o n gal axy formation models 
(e.g., lEisenstein fc Loebl j TO Ol iM cGaugh & de Blo3 
1998i iM o et al.' ^998"; "Cmirteau & Rix 199^ 
Steinmetz fc Navarro 1999; Courteau fc Rix 1999; 



van den BoschI l2000t iBullock et all l200lt iMo fc Mad 



20041 iMaverfcMoorel 120041: iGnedinetal 



Governato et al.l 120071: iTruiillo-Goniez et al 



de Rossi et al.l l2010t iTonini et al.l l2010t ). In a broader 



context, it provides a useful test of theories that 
seek to modify gravity in order to obviat e the need 
for dark matter (e.g iMilgromI Il983: Sanders! 119841: 
iMannheim fc Kazanaj Il989t lMoffatri2006f ). These 
should make specific and testable predictions for the 
Tully-Fisher relation. 

The Tully-Fisher relation was originally posed as an 
empirical relation between optical luminosity and the 
width of the 21 cm line. These are presumably prox- 
ies for more physical quantities. The 2 1 cm line-width 
is a proxy for rotation velocity (e.g., iVerheiienl 120011 : 
lYegorova fc Saluc"c]| 120071 ) which is widely presumed to 
be set by the dark matter halo. Luminosity is a proxy 
for stellar mass, which is presumably related to the to- 
tal mass of the system. Thus a relation like Tully-Fisher 
seems natural, with the mass of the parent halo specify- 
ing both luminosity and rotation velocity. 

Luminosity is not a perfect proxy for mass. The 
stellar mass-to-light ratio can, and presumably does, 
vary with galaxy type. This leads to Tully-Fisher re- 
lations with different slopes and scatter depending on 



the band-pass used to measure the luminosity (e.g. , 
TuUv et al.' 1998; Vcrheiien 200l iCourteau et all [20071: 
Noordermeer fc Verheiien 2007: [Masters et al.l l2008[ l. 
Moreover, low luminosity galaxies tend to fall be- 
low the extrapolation of the fit to bright galaxies 
(jPersic fc Saluccilll991al: [Matthews et al.lll998D . Yet ah 
these observed relations presumably stem from a single 
u nderlying physic al relation. 

[Freeman! (fT999( ) and [McGaugh et al.[ (|2000( ) suggested 
that baryonic mass is a more fundamental quantity 
than luminosity. For any given stellar population, 
different mass-to-light ratios will be produced in 
different bands by the same mass of stars. More 
importantly, physics should not distinguish between 
mass in stars and mass in other forms. Adding the 
observed gas mass to the stellar mass results in a 
Baryonic Tully-Fisher Relation (BTFR) that is linear 
(in lo g space) over m a ny dec a des in mass ( Vcrhciiefl 
tooT; 'Bell fc de Jong" '200T; 'Gurovic h etal.1 12004 



McGaugh 2005b; Pfcnniger fc Revaz 2005!: [Begum et al J 
2008at '[st ark et al.: [2009t [Trachternach et al.! I200'a 
Gurovich et al.. .2010l ). The BTFR appears to be the 



fundamental physical relation underpinning the empiri- 
cal Tully-Fisher relation. As such, the true form of the 
BTFR is of obvious importance. 

In this paper, I calibrate the BTFR with gas rich galax- 
ies. The masses of gas rich galaxies can be more accu- 
rately measured than those of the star dominated galax- 
ies that traditionally define the Tully-Fisher relation be- 
cause they are less affected by the systematic uncertainty 
in the stella r mass-to-light ratio jMcGaugh et al.l [2000 : 
[Stark et al.! [2009 : McGaugh 2011) . This results in a 
purely empirical calibration with a minimum of assump- 
tions. 

In §2, I discuss the data required construct the BTFR, 
assemble samples of gas rich galaxies, explore the statis- 
tics of the data, and calibrate the relation. In S3 I ex- 



amine the implications of the relation for ACDM and 
MOND. I define a feedback efficacy parameter £ to quan- 
tify the mapping between observations and theoretical 
dark halo parameters. I estimate the value of the MOND 
acceleration constant ag from the empirical BTFR. A 
summary of the conclusions are given in §4. 

2. THE BARYONIC TULLY-FISHER RELATION 

In order to construct the BTFR, it is necessary to esti- 
mate the baryonic mass of a galaxy on the one hand, and 
its rotation velocity on the other. The baryonic mass is 
the sum of all observed components, stars and gas: 



Mb = A/* + Ma 



(1) 



Combining this with measurements of the rotation veloc- 
ity, one obtains the BTFR: 



Mh^AV". 



(2) 



Remarkably, this simple formula suffices to describe 
the data for rotating galaxies. There is no indication 
of any need for a third pararneter such as radius or 
surface brightness (iZwaan et al.lll995t iSprayberrv et al.l 



^ ISprayf 
?hl20n5all : 



[ToMICcmrteau fc Rixlll999l : IMcGaughi20n5an : theTullv 
Fisher relation is the optimal projection of the fundamen- 
tal plane of disk galaxies. 

2.1. Stellar Mass 

The stellar mass of a galaxy is determined from its 
measured luminosity and an estimate of its stellar-mass- 
light ratio T,: 

M, = T^L. (3) 

The mass-to-light ratio is common ly obtaine d from 
stellar population models (e.g., iBell et all 120031 : 
iPortinari et al.l 120041 ) . These predict correlations be- 
tween T* and color so that the observed color is used to 
estimate T,. 

The use of a stellar mass-to-light ratio estimated from 
stellar population models causes a significant systern - 
atic uncertainty in the stellar mass (jConrov et al.ll2009l ). 
which is particularly sensitive to the IMF. Depending on 
the choice of model, one can derive BTFRs with slopes 
anywhere in the range of a; = 3 to 4 (McGaugh 2005 q). 
This is a basic limitation that is hard to avoid in sam- 
ples of spiral galaxies where most of the baryonic mass 
is in the form of stars. This can be overcome by use of 
samples of galaxies w hose baryonic mas s is dominated by 
gas rather than stars (jStark et al.ll20Q9l) . In such gas rich 
galaxies, the uncertainty in stellar mass is reduced to a 
minor contributor to the total error budget. Neverthe- 
less, I am careful to employ a self-consistent population 
model to all galaxies employed in the analysis (Table [1]). 

2.2. Gas Mass 

Galaxies with Mg > Af* can be found among very late 
type (Sd and Sm) spirals and Irregulars with rotation 
velocities Vf < 90 kms~^ (e.g.,|E der fc Schombert"2000': 
Walter et all 120081: IBegum et al. 2008b; Swatcrs et aL 



20091) . With enough such galaxies, it is possible to obtain 



an absolute calibration of the BTFR that i s effectively 
indep endent of the stellar mass estimator (jStark et al.l 
I2009D . Here I describe gas mass estimates for galaxies 
used in the analysis. 



2.2.1. Atomic Gas 

Neutral atomic (HI) gas usually dominates the gas 
component of disk galaxies. The mass of HI follows di- 
rectly from the physics of the spin-flip transition of hy- 
drogen, the cosmic hydrogen fraction, and the distance 
to each galaxy. There is no uncertain mass-to-light ratio 
as with stars. 

The gas mass is estimated by 



Mg = tjMhi, 



(4) 



where Mhi is the observed atomic gas mass. The factor 
77 accounts for the cosmic hydrogen fraction and forms of 
gas other than HI. Since the majority of objects of inter- 
est here are low luminosity late type galaxies, a primor- 
dial hydrogen fraction is assumed (X"^ — 1.33) rather 
than the commonly assumed solar value {X~^ = 1.4). 
This difference is less than the typical uncertainty in 
Mhi, but does make the numbers here differ slightly from 
those in the source data. 

2.2.2. Molecular Gas 

After atomic hydrogen, molecular gas is the next sig- 
nificant gas component. Molecular gas persists in be- 
ing difficult to detect in late type, low surface brightness 
galaxies. Where detected, it appear s to contain rela- 
tively l ittle of the int e rstell ar gas mass ("Matthews & Gai 
20011 [OTfciletal.1 [2OOOI: lO 'Ncil & Schinnerer 200" 



Matthews et al.ll2005i iDas et al.l I200a 12010|) . Such di 
rect measurements of molecular gas are not available for 
most galaxies. To estimate the molecular gas mass, we 
might resort to scaling relations for Mh2 /Mhi such that 



' X \ Mhi 



(5) 



lYoung fc Knezekl (|1989D showed that Mh^/Mhi depends 
on morphological type, albeit with large scatter. From 
their data, McGaugh fc de Blok (199 7) derived the scal- 
ing relation MhJMhi = 3.7 - 0.8T -I- 0.043T2. The 
vast majority of of gas rich galaxies are very late types 
(T > 8), so the typical correction is small (< 5%). 

We can u se the recent observations of THINGS galaxies 
([Walter et a l. 2008) to provide a n independent estimate 
of MhJMhi- iLerov et all (|2008f) use the observed (low) 
star formation in dwarf Irregular galaxies to infer the 
amount of molecular gas present. This maps smoothly 
into the observed molecular content of brighter galax- 
ies, with all disks following a relation between surface 
densities such that T,h2/^hi = S*/(81 Mopc"^) (their 
equation 33). This is a relation between local surface 
densities, while a relation between global masses is de- 
sired here. The atomic gas is typically spread over a 
much larger area than the stars, so we cannot simply as- 
sume that the ratio of surface densities is the same as 
the ratio of masses. However, the si ze of the stel l ar an d 
molecular disk are usually similar (jHelfer et al.l l2003| ) . 
Assuming that they are the same gives 



M, 



Ho 



^HI 



M. 



81M0PC- 



(6) 



In this case, the amount of molecular mass depends on 
the stellar mass rather than the atomic gas mass and 
morphological type as in equation ([5]). 



In order to apply equation ^, we need an estimate 
of the average HI surface density. This is almost al- 
ways < 10 MfTi pc~ ^ in low surface brightness galaxies 
(|de Bloket al.lll996( ). In the pr ototypical ca se DDO 154, 
T,Hi « 7 Mopc-2 at ah radii (jLerov et al.|[2008[ ). Tak- 
ing this value as typical implies that molecular gas is 
perhaps 9% of the stellar mass of these systems. For 
galaxies with less than half of their mass in the form of 
stars, we again infer that the molecular gas makes up 
< 5% of the baryonic mass. 

For either of these estimators, the molecular gas con- 
tent of Hl-dominated late type galaxies is small. So 
small, in fact, that it is less than the uncertainty in the 
measurement of the HI mass in most cases. This would 
seem to be a good working definition of negligible, so 
rather than make an uncertain correction for the molec- 
ular gas mass, we neglect it entirely. In effect, we assume 
a constant rj = 1.33. 

2.2.3. Other Forms of Gas 

We restrict our baryonic mass budget to the compo- 
nents that are actually detected, as discussed above. It 
is of course possible that other reservoirs of baryons exist 
in some undetected form. There is, however, no direct 
evidence for baryons in any form that is comparable in 
mass to the stars and atomic gas. 

One conceivable reservoir of baryonic mass is hot gas in 
the halo. Such a component is expected in some galaxy 
formation scenarios. While some hot gas is certainly 
present, it appears that its densi ty is too low to make 
a sub stantial contribution here (jAnderson &: BregmanI 

IMol) . 

Another suggestion is that some of the dark matter 
could be in t he form of undetected c old molecular gas in 
galaxy disks (jPfenniger et al.l[l994( ). iPfenniger fc Reva3 
HoOJI) showed that a marginal improvement i n the scat- 
ter of the BTFR of iMcGaugh et~al] pOOOt ) could be 
obtained by treating rj as a, free parameter, and sug- 
gested that the best fit value (?/ w 3) implied addi- 
tional dark baryonic mass in the disk. In more accu- 
rate data, the scatter is consistent wi th being caused 
entirely by observational uncertainties fVerheijen' '2001t 
McGaugh 2005b; Sta rk ct al. 2009; Trachtcrnach ct alj 
2009HMcGaughll20lH) . nullifying this motivation for dark 



baryons. 

In sum, there is no compelling evidence for substantial 
reservoirs of baryons in disk galaxies in any form other 
than those that are directly detected. I therefore con- 
sider only the physically motivated value of ?7 = 1.33. 
Obviously, a different result would be obtained if we al- 
lowed v to include hypothetical mass comp onents (cf. 
IPfenniger fc Revazj[2005L iBegum et al.ll2008al ). 

2.3. Rotation Velocity 

The characteristic circular velocity of a disk galaxy can 
be measured with resolved rotation curves, or approxi- 
mated by line- widths. The latter is a flux- weighted inte- 
gral over the velocity field, and lacks the clear physical 
meaning of a rotation curve. It is, however, readily ob- 
tained from single disk 21 cm observations, whereas ob- 
taining and analyzing interferometri c data cubes is con- 
siderably more effort intensive fe.g.. lde Blok et al]|2008l : 
pKachternach et al.„2009,V 



When a full rotation curve is available, one can make 
different choices about where to measure the circular 
velocity. Common choices i nclude the maximum ob- 
serve d velocity (Vmax, e.g., iNoordermeer fc Verheijenl 
l2007f l. that in the outermost measured regions where 
rotation curves tend towards b e coming approxima tely 
flat {Vf, e.g., IVerheiienI [200ll : IMcGaughl l2005bD or 
the velocity at some p a rticular optical radius (e.g. , 
Persic fc Saliicc| Il99m IMcGaugh fc de Blokl IT998al: 
Courteau fc Rixlll999( ). The gas rich galaxies of great- 



est interest here typically have gradually rising rotation 
curves that slowly bend over toward becoming flat with- 
out ever having an intermediate peak. For this com- 
mon morphology, all these measures of velocity are very 
nearly equivalent. In particular, Vmax ~ Vfi and the ve- 
locity measured at a particular optical radius is also the 
same unless that radius is chosen to be very small (see 
r^gorova fc Salucci 2007 ). Indeed, it is not particularly 
meaningful to define an optical radius for gas rich galax- 
ies, as the scale length of the baryonic mass c an be con- 
siderably larger than that of the s tars alone (jMcGaughl 
I2005at ISeii^ood fc McGa^Ighll2005h . 

For quantitative analysis, I restrict consideration to 
cases where Vf is explicitly measured, simply adopting 
the measured velocity from each source. For compari- 
son purposes, I also show line-width {W20) data, but do 
not include these data in the calibration of the BTFR as 
line-widths and rotation velocities are comparable but 
not identical measures. In order to make comparisons, 
I assume Vf = W20/2 for gas dominated galaxies with 
line-width measurements on the presumption that they 
would show rising rotation curves like those observed for 
comparable galaxies. For the star dominated galaxies, 
I assume W20/2 = Vp = 1-lV/. This qualitatively ac- 
counts for the fact that bright galaxies have rotation 
curves that peak before falling gradually toward flatness 
(jPersic fc Salucci|[l991bD . The quantitat ive factor is cho- 
sen to reconcile the BTFR fit to W20/2 (jMcGaugh et al.l 
[loM, ^ = 35 Mq km"^ s*) with that fit to Vf (jMcGaughl 
12005b, A = 50 Mq km~^ s^*). This is only done to facili- 
tate comparison between independent samples. 

To account for non-circular motions, I accept what- 
ever prescription was used by each source. In gen- 
eral, these correction s turn out to be rather small (cf. 
TuUv fc Fouaue 1985tlOh et al.[[200llTrachternach et al.l 






2009t iDalcanton fc Stild I2010D . as the typical veloc- 
ity dispersion of lo w mass galaxies is ~ 8 kms""'^ 
(jKuzio de Narav et al. 2009). Since the correction to the 
circular velocity occurs in quadrature, it is only some- 
what important only in the slowest rotators. There ex- 
ist dynamically cold, rotating disks down to at least 
Vf ^ 30 kms~^, and per haps as low as ^ 10 kms~^ 
(jTikhonov fc Klvninl[2009h . 

2.4. Gas Rich Galaxy Samples 

There is an enormous amount of data in the literature, 
but only a finite amount of it pertains to gas dominated 
galaxies with Mg > Af * that also have all the required in- 
formation for constructing the BTFR. I describe here the 
accumulated data to date. I also include associated sam- 
ples of star dominated galaxies. This is done strictly for 
comparison as the uncertainties for star dominated galax- 
ies are dominated by the choice of stellar IMF. This is 



the systematic uncertainty that we ehminate by making 
use of gas dominated galaxies. For gas rich galaxies, the 
major sources of uncertainty are the 21 cm flux measure- 
ment and the distance to each galaxy. The uncertainty 
in the measurement of Vf is usually much smaller than 
that in mass (Table [Ij , provided that the inclination is 
reasonably well constrained. Cases where this condition 
might not be satisfied are discussed and, if suspect, ex- 
cluded from fitting. 

2.4.1. Stellar Mass Estimates 



T he data oflBegum et all (l2008aD . iStark et all (|2009t ). 
and 'Trac hternach et al.l ( 2OO90 constitute the sample of 
gas rich galaxies for which all the necessary information is 
available to construct the BTFR. For these galaxies, stel- 
lar mass is, by definition, a minority of the baryonic mass. 
Nevertheless, we take care to estimate stellar masses self- 
consistently. For specificity, we use their observed B — V 
colors and the populati on mode l of Portinari c t al.l ()2004f ) 
with a Kroupa IMF ()Kroupal [2002[ ). Other plausible 
choices make no difference to the r esult because the gas 
mass dominates (jStark et al.ll2009[ ). 

Several other samples are inclu ded for comparison 
(iMcGaugh et al.l 12000: McGaugh 200lB iGurovich et"all 
l201Clf ). Since these data are not used in the analysis, we 
make no effort to reconcile them to the same population 
model and simply adopt the stellar mass estimated by 
each source. Despite the minor differences in the meth- 
ods used to estimate stellar mass, these data are remark- 
ably consistent with one another. 

2.4.2. Galaxies with Rotation Curve Data 



\Beaum et al\ i2008d) : — Data for remarkably slowly 
rotating gala,xies a re presented and analyzed by 
iBegum et al.l (|2008ai rbl). The corrections for non-circular 
motions are substantial in some cases. However, the 
most significant difference in the analysis here is that 
we determine the gas mass as discussed in i)2.21 and 
we select galaxies for consistency in their inclination 
determinations. It is not uncommon for the opti- 
cal portions of gas rich galaxies t o be much sinaller 
than their dominan t HI disks (e.g., iMeurer et all 119961 : 
iGentile et al.l I2007D and for the outer HI isophotes 
to have a different shape than the opt ical isophotes 
(|0h et all l200llTrachternach et al.l 120091 ). If the stars 
suffer an oval distortion (bar), then an inclination based 
on optical isophotes will be systematically over estimated 
(IMihos et a l.lll997l: |de Blok fc McGaughTfTOOl . Indeed, 
IBegum et al. (2008^ observe that the optically inferred 
inclinations often exceed those inferred from HI. One 
must be wary that beam smearing might make the HI 
isophotes rounder than they should be, potentially caus- 
ing an underestimate of inclination. To avoid either 
effect, I demand consistency between optical and HI 
inclination determinations. Specifically, I require that 
sin(iopt) be within 12% of sin(z///). The restriction is 
placed on sin(z) because this is what matters to the ac- 
curacy with which Vf can be measured. We are con- 
cerned with the velo city at large radii, so the HI inclina- 
tion is adopted (as in IBegum et al.ll2008af) : the difference 
iHi I is treated as an additional source of er- 



the need to retain at least some objects. W ith this crite- 
rion, 16 objects are retained from Table 1 of lBegum et al.l 
(|2008al) and 12 are rejected. Requiring any greater ac- 
curacy eliminates nearly the entire sample. UGC 5456 
is excluded for being marginally st ar dominated. T hree 
objects are also in the sample of iStark et al.l ()2009[ ). 



\ Stark et all 1200 A ): — Data for gas rich galaxies are se- 
lect ed from the l i teratu re with an emphasis on data qual- 
ity. IStark et all (|2009[ ) required resolved rotation curves 
that were observed to flatten sufficiently that the log- 
arithmic slo pe became < 0.1 in their outer parts (see 
iSwaters et al. 2009. from which much of this sample orig- 
inates). This requirement prevents the data from skew- 
ing to low Vf due to the inclusion of galaxies whose ro- 
tation curves are still rising rapidly at the last measured 
point. 



I Trachternach et al\ (MM):— New WSRT observations of 

II very gas rich galax ies selected from the sam ple o f 
iSchombert eTall (|199l and lEder fc Schombertl lJ2000l ) 
are used to measure both W20 and Vf. Of the 11 galax- 
ies observed, 6 show clear evidence of rotation. Interest- 
ingly, these all have HI profiles that are double-horned 
or at least show some structure. Galaxies for which rota- 
tion is not positivel y detecte d tend to have single-horned 
profiles. Trachternach et al.l (i2009.) discuss two methods 
for measuring the outer circular velocity that are largely 
consistent. I adopt for Vf their values measured from 
the position- velocity diagram (Vpy from their Table 3). 
For D575-2, I estimate V/ = 59 ± 7 kms"^ based on the 
outermost detected 21 cm emission in their data. 



\McGauoli l[2005U) : — The data for star dominated galax- 
ies with Vf measured from extended HI rotation curves 
are shown for comparison. They are not used in the 
analysis here. For illustration, stellar masses adopt the 
Q = 1 MOND mass-to-light ratios. The BTFR fit to 
these data [Mb = (50 M© kni'^s")!^/] can be tested by 
extrapolation to the much lower circular velocities of the 
new gas rich galaxies. 

The thr ee data sets of IBegum et all (l2008af) . 

IStark et all (|2009l ).and [Trachternach et al.l(|2009l ) are the 
most interesting for our purposes here. They compose a 
sample of gas dominated galaxies with resolved rotation 
curve measurements. These data are listed in Table [1] 
and arc used in the analysis. The data for these galaxies 
is less precise than those for bright star dominated galax- 
ies. However, they are more accurate, as their masses are 
not sensitive to systematic uncertainties in the IMF and 
stellar population mass-to-light ratios. 

2.4.3. Galaxies with Line-width Data 

Samples of galaxies with line-width data are used for 
comparison purposes only. We shall see that line-width 
data give a systematically shallower slope and greater 
scatter for the BTFR than do rotation curve data. This 
appears to be a limitation of the line-width as a proxy 
for the rotation velocity. Aside from this detail, the data 
present a consistent overall picture. 



Ai 



^opt 



ror. The particular threshold of 12% is determined by 
trial and error, balancing the desire for accuracy with 



\McGauah et al.\ l(2000i) : — Gas rich galaxies wer e se- 
lected fro m the sample o f ISchombert et al.l (|1997f) and 
lEder fc Schombertl (|2000f ) to have inclinations i > 45°. 
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Figure 1. Baryonic mass as a function of line- width (left) and circular velocity as measured from resolved rotation curves (right). Left 
panel: data from McGaugh ct al. (2000) are shown as light gray points with uncertainties suppressed for clarity. Gas dominated galaxies 
(with M* < Mg) are shown as small circles and star dominated galaxies (Af* > Mg) are shown as small squares. More recent data are 
shown with larger symbols. The data of Trachternach ct al. (2009) are shown as dark blue stars in cases where rotation is detected and 
as open stars when it is not. The data of Begum et al. (2008a) are shown as light blue squares for cases with consistent optical and HI 
inclinations and as crosses in cases where these differ substantially. The H I data oflGurovich ct al. (2010) are shown as orange left-pointing 
triangles. The data of the star dom inated galaxies of Sakai et ah] |2000l ) are shown as red right-pointing triangles with stellar mass from 
the H-hand data of Guro vich et all ([2010). The two samples of star dominated galaxies are compared in the inset. This also serves to 
illustrate the range covered by the bright spirals that traditionally define the TuUy-Fisher relation — a range that is vastly expanded by 
the inclusion of gas rich dwarfs. In the right panel, data for galaxies with Vf measured from resolved rotation curves include the rotating 
cases of Tra chternach et al. (2009), the data of Begum et al. (2008a) with con sistent incl i nations , the gas dominated galaxies compiled by 
IStark et al.l l|2009l . green circles), and the star dominated galaxies compiled bv lMcGaughl H2005bl . dark gray squares). 



Stellar masses were estimated from /-band data assmii- 
ing a constant mass-to-light ratio (T^ = 1.7 Mq/Lq). 
We retain this original prescription, making on ly the 
correction to stellar mass discussed by Gurovich et al.l 
([2004). Data for star domina ted galaxies w i th JJ- band 
photometry were taken from iBothun et al.l ()1985D . A 
constant H-hand mass-to-light ratio of T^ = 1 Mq/Lq 
was assumed. These data are included here for compari- 
son to more recent data but are not used in the analysis. 



New data are presented for 
Detailed stellar population 



\Gurovich eTai} \201(\) : — 
galaxies with Mg > M^ 
modeling has be en performed to dete rmine the stellar 
mass (Table 11 of lGurovich et al.ll2010[ ). The same mod- 
eling procedure has also been applied to the star domi- 
nated sample of Sakai ct al. (2000). We adopt the stellar 
mass determinations based on iJ-ba nd photom etr>{3 for 
these galaxies (Table 4 of lGurovich et al.ll201Q) to facil- 
itate direct comparis on with the data of Bothun et all 
(fT985l) employed by (McGaugh et alJ|2000D. AH of the 
bright galaxies of Sakai ct al. (2000) have nicely double 

^ Whether optical or JJ-ba nd data are used to estimate the stellar 
mass for the bright galaxies of I Sakai et al.l I I2000I) makes a difference 
to the slope of the BTFR: a; = 3 in the optical and a; = 4 in 
the //-band. This exemplifies the systematic uncertainty that gas 
dominated galaxies save us from. 



horned line profiles with well measured line- widths. This 
is not the case for the n ew gas rich galaxies presented 
by IGurovich et al.l (|2010) , for which the HI profiles are 
all narrow and single-horned (Gurovich, private commu- 
nication). Such profiles can be poor indicators of Vf 
(|Trachternach et al.ll2009l ). 

2.5. Comparison of the Data 

Fig.|T]shows the data. Over all, the agreement between 
the various data is good. This is true even for the star 
dominated galaxies, despite the differe nces in the popula - 
tion modeling. The spiral gala xies oflSakai et al.l (I2OOO0 
fall in the midst of the data of IBothun et al.l ( 19851 ) fin- 
set of Fig. [Ij with indistinguishable slope and normaliza- 
ti on. The chie f difference is the lower scatter in the data 
of lSakai et al.l ([2000,). This presumably results from the 
more accurate distanc es to the individual ga laxies based 
on HST observations (iFreedman et al.ll200lD . This is ex- 
pected: as the data improve, the scatter is reduced. In- 
deed, essentially all of the scatter can be accounted for 
by observational uncertainty and the expected variation 
in stel l ar mass-to-light r atios (|Verheiienll2001l : iMcGaughl 
l2005btlStaFk et al.ll2009l) . 

The scatter is larger for lower luminosity galaxies. This 
is also expected: it is more challenging observationally to 
obtain precise data for faint dwarfs than it is for bright 
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Figure 2. The data in a coordinate system rotated (equations [7] and (Sj so that the abscissa is the position of a galaxy along the BTFR 
and the ord inate is its perpendicular distance from the BTFR. Symbols as per Fig.[T]with the exception that the star dominated galaxies of 
ISakai et al.l IJ2000) are shown t wice: once with the H-band stellar masses of Gurovich et al. (2010) as in Fig. [l] and again for the hybrid V 
and H band stellar masses of lGurovich et al] 1 12010 . open right-pointing triangles). This illustrates the systematic uncertainty that plagues 
star dominated galaxies but from which gas dominated galaxies are largely immune. The top panel shows the BTFR with slope 3 fit to the 
data of (Gurovich et al.ii 20101 '). The bottom panel shows the BTFR fit of lMcGaughl (,2005b) with slope 4. The data shown in the bottom 
panel are independent of those to which the fit was made. 

spirals. Nevertheless, the data for the gas rich dwarfs 
appear to be largely consistent. Moreover, the scatter 
is reduced as we progress from line-widths to resolved 
measures of the rotation velocity (left vs. right panels in 
Fig.[T]). The improved accuracy with which the rotation 
velocity is measured is reflected in reduced scatter. 

There are some outliers among the gas rich dwarfs, 
especially when W20 is used as the velocity estimator. 
The two most dram atic outliers from the sample of 
IBeeum et al.] (l2008al) are UGC 3755 and KK98-65, which 
fall far to the left of the main correlation. These are re- 
jected by the requirement of consistent optical and HI 
inclinations as are a number of galaxies that do lie along 
the BTFR in Fig. [T] when the HI inclina tion is used. The 
outliers from lTrachternach et al.l (|2009l ) are those objects 
for which rotation is not clearly detected. (D575-5 is 
the m ost discrepant c a se). T hree objects from the sam- 
ple of lGurovich et all ()2010D stand out from the rest of 
the data: ESQ 085-G 088, AM 0433-654, and HIPASS 
J1424-16b. ESQ 085-G 088 has a large velocity error, 
so its status as an outlier is only marginal. The other 
two objects clearly deviate from the rest of the data, in - 
cluding the bulk of the data of iGurovich et al.l ()2010[ ). 
Since we do not have resolved measures of rotation for 
these objects, it is possib le that they, hke the single- 
peaked HI profiles of Trachternach et al.l ()2009f ). would 
not show clear evidence of rotation if so observed. An- 
other possibility is th at the optical inclinations used by 
IGurovich et al.l ()2010f) are systematically in error as dis- 
cussed in i)2.4.2l Lacking resolved HI observations, we 
can not make the same consistency check that is possible 



with the data of iBegum et al.l (j2008a). The inclination 
of AM 0433-654 is 43° and that of HIPASS J1424-16b is 
47°. A shift of ^ 0.15 dex in VF20 would reconcile these 
galaxies with the other data. If inclination is responsi- 
ble, the corresponding inclinations would be 29° for AM 
0433-654 and 31° for HIPAS S J1424-16b. Fro m cxpcri- 
ence with such objects ( e.g., Ide Blok fc McGa ugh 199^; 
'Trachtern ach et ani2009f l. this is well within the realm of 
possibility. Indeed, it would be surprising if something 
like this did not happen with only optical inclinations to 
work with. 

Over all, the data paint a consistent picture (Fig. [ij 
despite the inevitable concern over the accuracy of a few 
individual measurements. The BTFR is continuous and 
linear in log space over five decades in mass. There is a 
hint of curvature when W20 is used as the rotation veloc- 
ity measure, but it is not statistically significant. There 
is not even a hint of curvature when resolved rotation 
curves are available to measure Vf. 

We expect the accuracy of the data to improve as we 
move from integral measures like the line- width to phys- 
ically meaningful rotation velocities measured from re- 
solved HI data cubes. This expectation is realized in the 
reduced scatter seen in the right panel of Fig. [TJ We can 
quantify this with ordinary least-square fits to the data. 

If we app l y the maximum likelihood method of 
iWeiner et al.l (|2006f ) to fit the new line- width data (i.e., 
all the data in the le ft panel of Fig. [1] except those of 
lMcGaugheraIll2000D . we obtain a slope x = 3.41 ± 0.08 
and an intercept log^ = 2.91 ± 0.17. If we perform 
the same exercise with the rotation curve data in the 



Table 1 

Data for Selected Gas Rich Galaxies 



Galaxy 


D 


Vf 


log(M.) 


log(Mg) 


Ref. 




Mpc 


kms ^ 


Mq 


M© 




DDO 210 


0.94 


17 ±4 


5.88 ±0.15 


6.64 ± 0.20 


1,2 


Cam B 


3.34 


20 ±12 


6.99 ±0.15 


7.33 ±0.20 


1,2 


UGC 8215 


4.5 


20 ±6 


6.81 ±0.15 


7.45 ±0.20 


1 


DDO 183 


3.24 


25 ±3 


7.24 ±0.15 


7.54 ±0.20 


1 


UGC 8833 


3.2 


27 ±4 


6.94 ±0.15 


7.30 ±0.20 


1 


D564-8 


6.5 


29 ±5 


6.76 ±0.20 


7.32 ±0.13 


3 


DDO 181 


3.1 


30 ±6 


7.26 ±0.15 


7.56 ±0.20 


1 


P51659 


3.6 


31 ±4 


6.67 ±0.15 


7.85 ±0.20 


1 


KK98 246 


7.83 


35 ±6 


7.72 ±0.15 


7.93 ±0.20 


1 


UGCA 92 


3.01 


37 ±4 


7.78 ±0.15 


8.32 ±0.20 


1 


D512-2 


14.1 


37 ±7 


7.58 ±0.20 


7.96 ± 0.06 


3 


UGCA 444 


0.95 


38 ±5 


7.34 ±0.15 


7.75 ±0.20 


2 


KK98 251 


5.6 


38 ±5 


7.34 ±0.15 


8.02 ±0.20 


1 


UGC 7242 


5.4 


40 ±4 


7.57 ±0.15 


7.78 ±0.20 


1 


UGC 6145 


7.4 


41 ±4 


7.20 ±0.15 


7.56 ±0.20 


1 


NGC 3741 


3.0 


44 ±3 


7.24 ±0.15 


8.45 ±0.20 


1,2 


D500-3 


18.5 


45 ±6 


6.97 ±0.20 


7.94 ±0.05 


3 


D631-7 


5.5 


53 ±5 


6.88 ±0.20 


8.29 ±0.15 


3 


DDO 168 


4.3 


54 ±3 


8.07 ±0.15 


8.74 ± 0.20 


2 


KKH 11 


3.0 


56 ±5 


7.28 ±0.15 


7.85 ±0.20 


1 


UGC 8550 


5.1 


58 ±3 


8.25 ±0.37 


8.46 ±0.39 


2 


D575-2 


12.2 


59 ±7 


7.63 ±0.20 


8.62 ±0.07 


3 


UGC 4115 


7.5 


59 ±6 


7.77 ±0.15 


8.58 ±0.20 


1 


UGC 3851 


3.2 


60 ±5 


8.45 ±0.15 


9.09 ±0.20 


2 


UGC 9211 


12.6 


64 ±5 


8.12 ±0.39 


9.21 ±0.41 


2 


NGC 3109 


1.3 


66 ±3 


7.41 ±0.15 


8.79 ±0.20 


2 


UGC 8055 


17.4 


66 ±7 


8.09 ±0.15 


9.02 ±0.20 


1 


D500-2 


17.9 


68 ±7 


7.41 ±0.20 


9.06 ±0.05 


3 


IC 2574 


4.0 


68 ±5 


8.94 ±0.15 


9.20 ±0.20 


2 


UGC 6818 


18.6 


72 ±6 


9.22 ±0.16 


9.28 ±0.20 


2 


UGC 4499 


13.0 


74 ±3 


8.75 ±0.27 


9.32 ±0.29 


2 


NGC 1560 


3.45 


77 ±3 


8.70 ±0.15 


9.23 ±0.20 


2 


UGC 8490 


4.65 


78 ±3 


8.36 ±0.15 


8.96 ±0.20 


2 


UGC 5721 


6.5 


79 ±3 


8.17 ±0.37 


9.05 ±0.39 


2 


F565-V2 


48. 


83 ±8 


8.30 ±0.21 


9.04 ±0.24 


2 


F571-V1 


79. 


83 ±5 


9.00 ±0.19 


9.33 ±0.22 


2 


IC 2233 


10.4 


84 ±5 


8.96 ±0.15 


9.32 ±0.20 


2 


NGC 2915 


3.78 


84 ±10 


7.99 ±0.15 


8.78 ±0.20 


2 


NGC 5585 


5.7 


90 ±3 


8.98 ±0.38 


9.27 ±0.40 


2 


UGC 3711 


7.9 


95 ±3 


8.92 ±0.15 


9.01 ±0.17 


2 


UGC 6983 


18.6 


108 ±3 


9.53 ±0.16 


9.74 ±0.20 


2 


F563-V2 


61. 


111 ±5 


9.41 ±0.17 


9.63 ±0.21 


2 


F568-1 


85. 


118 ±4 


9.50 ±0.18 


9.87 ±0.22 


2 


F568-3 


77. 


120 ±6 


9.62 ±0.18 


9.71 ±0.22 


2 


F568-V1 


80. 


124 ±5 


9.38 ±0.18 


9.65 ±0.22 


2 


NGC 2403 


3.18 


134 ±3 


9.61 ±0.15 


9.77 ±0.20 


2 


NGC 3198 


14.5 


149 ±3 


10.12 ±0.15 


10.29 ±0.20 


2 




iS. — 1 










Reference 


IBeeum et al.l Il2008all: 2. 


IStark et al.l ^200fll): 


3. ITrachternach et al 


1 J2OO9I). 








Note. — 


- An electronic version of this table is available at 


|http://www. 


astro.umd.edu/"ssm/data/gasrichc 


atatable.txt 1 





right panel of Fig. [I] we obtain x = 3.94 ± 0.11 and 
log A = 1.81 ± 0.22. These formal fits have rather dif- 
ferent slopes even tho ugh the chief differe nce is in the 
scatter. The method of lWeiner et al.l (|2006l ) accounts for 
uncertainties in both velocity and mass, and can option- 
ally estimate the amount of intrinsic scatter in a sam- 
ple. Accounting fo r the stated observational errors (e.g., 
iWeiner et al.ll200"6l ) implies an intrinsic scatter in the line- 
width data of ^ 0.2 dex. That is, there is 0.2 dex of 
scatter not explained by observational uncertainty. In 
contrast, this number is zero for the rotation curve data 
in the right hand panel: observational errors suffice to 
explain all the scatter, leaving no room for intrinsic scat- 
ter. This difference in scatter implies some source of 
variance in line-width measures that is perhaps not ac- 



counted for in the stated uncertainties, and not likely to 
be meaningful to the physics we are interested in here. 

At least part of the reason for the difference in fit- 
ted slopes is physical: W20 is not a perfect proxy for 
Vf . Bright galaxies typically have peaked rotation curves 
with Vmax > Vf. Line- width measurements do not re- 
solve where emission originates, and are thus sensitive to 
the maximum velocity rather than that at larger radius. 
This immediately skews the slope to shallower values. 
A similar if less pronounced effect can be seen in opti- 
cal rotatio n curves that are resolved but not particularly 
extended (jTorres-Flores et al.|[201lD . 

For the physical interpretation of the BTFR, we wish 
to have as precise a probe as possible of the quantity of 
interest. For ACDM, this would be the virial velocity of 
the dark matter halo. In the case of MOND, we need 
a measure that probes the low acceleration regime. In 
both cases, Vf provides a better measure than W2o- 

2.6. Statistics of the Gas Rich Galaxy Data 
2.6.1. Slope 

The slope of the BTFR has been addressed by nian:v 
authors fe. g.jMcGaugh e t al."2000^'Bell fc de Jon dl200li: 



Gurovich et al.l 120041: iM cGaugh 2005b; Begum et al.l 
2008ai IStark et al.l l2009t IGurovich et al.l 12010 ). and is 
typically found to be between x = 3 and 4. Most of this 
variation is due to the presc ription adopted fo r estimat- 
ing the mass-to- light ratio (jMcGaughl I2005bl ) . which is 
by far the strongest systematic effect. This can be reme- 
died by selecting gas dominated galaxies as described in 

^2xn 

IStark et all (|2009l ) fou nd a slope x = 3.94 ± 0.07 
with a bisector fit while IGurovich et al.l (J2010D found 
X = 3.0 ± 0.2 from their gas dominated galaxies (their 
Table 12). These result s appear to be contradi ctory. Re- 
peating the analysis of IGurovich et aLJ (|2010D . I obtain 
the same result. However, simply looking at the data, 
it is not obvious that they are in conflict. Part of the 
difference is simply the difference in velocity estimators, 
as line-widths give shallower slopes than resolved rota- 
tion velocities. Their sample is also small, so is easily 
influenced by a few outliers as discussed above ( ^2.5p . 

This difference between a; sa 3 and 4 is precisely the 
difference between the nominally anticipated slopes in 
ACDM and MOND. The correct answer is therefore of 
considerable importance and lends itself to controversy. 
Bear in mind that the data themselves are quite con- 
sistent (Fig. |T]), so any controversy is in the fitting and 
choice of velocity estimator, not in the data. The line- 
width is a convenient proxy for the rotation velocity when 
resolved rotation curves are unavailable, but not a substi- 
tute for them when predictive theories are being tested. 

To illustrate the difference between slopes of 3 and 4, 
we rotate the BTFR into a new coordinate system {i,S) 
parallel to it such that £ measures the position of a galaxy 
along the BTFR and S its deviation perpendicular to it. 
That is, 

^ = cos 61 log Vf + sin 6'(log Ah ~ log A) (7) 

(5 = cos 6i(log A/fc - log A) - sin d log Vf (8) 

where tan 6* is the slope of the BTFR. In Fig. [2] we 
show t he data rotated into the frames defined by the 
fits of .Gurovich et al.l (|2010. ) with slope tan6' — 3 and 
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Figure 3. Fits to tiie BTFR of gas dominated galaxies. Symbols as per Fig. [T] The solid lines show the best fit to the whole sample 

i with slope X = 3.82) and the bisect or (slope 3.98) of the forward and reverse fits (dashed lines) l imited to the samp les of Stark ct al] 
20091 ) and ITrachternach et al.l 1)20091 ) — that is, excluding the less accurate points (squares) from iBegum et al] l)2008ai ). The various fits 
are consistent within the uncertainties, and yield the same basic result. 

tively, with uncertainty ±0.24. The slight differenee that 
occurs when also fitting the data of Begum et al. (2008g) 
may result from a slight skew to low velocities, and tes- 
tifies to the difficulty i n obtaining quality data f or very 
slow rotators (see also ITrachternach et al.l [2OO90 . The 
skew effect becomes more pronounced if we relax the se- 
lection criterion that requires consistent optical and HI 
inclination determinations. Irrespective of these details, 
these fits are all consistent with one another (Fig. [J]), and 
none yield a slope that is meaningfully different from 4. 



iMcGaughl (|2005bf ) with tan6l = 4. 

A slope of 4 does a good job of describing the bulk 
of the data: there is no systematic deviation from the 
dotted line in the lower panel of Fig. [5] In contrast, there 
is an obvious residual from a slope of 3 in the the top 
panel. The data show a clear preference for the steeper 
slope. 

Fig. [5] also illustrates the systematic caused by dif- 
ferent mass-to-light ratio estimates for star dominated 
galaxies. The galaxies from lSakai et al.l (|2Q0Clf ) are shown 
tw ice, once with masses based on H-hand data (Table 4 
of iGurovich et al.l 120100 and again with the average of 
V an d H band stellar masses (Table 5 of IGurovic h et all 
I2010f) . The hybrid V-H masses are nicely consistent with 
a slope of 3, while H alone prefers a slope of 4. While one 
must choose which band to trust for bright galaxies, this 
choice is rendered irrelevant for gas dominated galaxies. 
We therefore concentrate the rest of our analysis on the 
gas rich galaxy sample described in S §2.4.2l and l2.4.ll 

Taking the 47 gas rich galaxies with reasonably trust- 
worthy data (TableU]), the best fit slope is x = 3.82±0.22 
with intercept log A — 2.01 ± 0.41. Various meth- 
ods (forward, reverse, and maximum likelihood) all give 
the sam e result. Of these da ta, the velocity measure- 
ments of IBegum et al.l (|2008aD are the most challenging: 
only three o f thei r objects meet the quality crit eria of 
I Stark et"aI1 ()2009l ). Excluding the data of Begu m~et al.l 
(I2008af) but retain i ng th ose of .Stark et al- (.20091 ) and 
ITrachternach et all ()2009D gives x = 4.05 ± 0.29 from a 
forward fit, 3.92 zfc 0.29 from a reverse fit, and a bisector 
slope (Iso be et al.| [l990') of 3.98 ±0.06. The correspond- 
ing intercepts are log^ = 1.57, 1.82, and 1.71, respec- 



2.6.2. Normalization 

An interesting question we can pose to the gas rich 
galaxy data is what the normalization of the BTFR is 
with slope fixed to 4. The gas rich galaxies provide an 
absolute calibration of the BTFR that is very nearly in- 
dependent of the estimator we use for the stellar mass- 
to-light ratio. By fixing the slope, we can sharpen our 
estimate of the intercept, and also examine higher order 
moments of the distribution. In addition to the scatter 
perpendicular to the BTFR, as , we can also examine the 
skew as and kurtosis a^. 

Table [5] shows these statistics for several combinations 
of the rotation curve samples. We start with the sample 
of lStark et al.l ()2009[ ) as their sample was constructed for 
this ex ercise. Adding the galaxies of 'Trachte rnach et al.l 
(|200 ^ decrease s the intercept slightly. Adding the data 
of B egum et al. (2008a) r aises a, gain, back near to the 
value found bv iMcGauehl (I2005bl . log A = 1.70). 

Table [5] gives quantitative estimates for the intercept. 
For each combination of samples both a weighted and un- 
weighted estimate is given. The weighted estimate comes 



from minimizing xt by adjusting log A with fixed slope: 
we minimize the orthogonal deviations. The unweighted 
estimate is made by finding the centroid of the distribu- 
tion in the orthogonal deviations 6 of the rotated data 
and requiring the centroid to lie at 5 = 0. The latter 
method gives no heed to the accuracy of the data, but is 
less sensitive to the influence of outliers whose formal 
uncertainty may understate systematic effects. These 
methods give indistinguishable results. 

All combinations of the various samples give consistent 
estimates for the intercept. To be specific, we adopt the 
weighted estimate from the combination of the Stark- 
Trachternach sample: 

A = 47±6MQkm"'^s'*. (9) 

This is the midpoint of the various determinations and 
has the smallest x^ ■ This value is of course very similar to 
the unweigh ted value (45 M0 km~ s^) obtained from the 
Stark et al! (2009) sample alone ()McGaugh et all 120101: 
McGaugh fc Woli l2QTQh . It is also similar to the value 
(50 Mq km"'' s'') found bv lMcGaughl (|2005ht ) from star- 
dominated galaxies, but has the virtue of being com- 
pletely empirical with little sensitivity to how stellar 
masses are estimated. 

It appears that at this juncture, systematic effects and 
the necessary assumptions are as limiting as data quality. 
It is traditional to assume rj — lA which is appropriate 
for solar metallicity atomic gas, but is slightly larger than 
the primordial 1.33 taken here for these little evolved, 
gas rich, low mass galaxies. The choice of rj = 1.4 is not 
importa nt to the star dominated galaxies in iMcGaugbl 
(|2005bD , but we note that adopting it here would cause a 
change in mass of at most 5% (for purely gaseous galax- 
ies) that would drive A = 47 ^- 49 M© km^** s**. This is 
at the level of our best guess in the mass of molecular 
mass. At this level the systematic uncertainty in the IMF 
that plagues stellar mass estimates also begins to play a 
role even in gas dominated galaxies. Given these con- 
siderations, any value in the range 45 — 50 Mq km^** s** 
is plausible, with values much outside this range quickly 
becoming less plau sible. If we u se the BTFR to estimate 
stellar masses ( Star k et al.ll2009l ). then more massive star 
dominated galaxies are consistent with a Kroupa IMF. 
Grossly different stellar masses would imply a break in 
the slope of the relation around ^ 100 kms~^ where the 
galaxy population transitions from frequent to infrequent 
gas domination. 

2.6.3. Scatter 

The scatter in our rotated coordinates is as = 0.06 
with only the third digit varying from sample to sample. 
In terms of mass this is 0.24 dex. The bulk of this uncer- 
tainty (0.20 dex) is provided by the error in the baryonic 
mass. Velocity uncertainties are less important in most 
samples, and only rival the uncertainty in mass in the 
lowest mass systems. Together, the uncertainties in mass 
and rotation velocity suffice to account for essentially all 
of the observed scatter. 

A robust estimator of the dispersion of a sample is 
the median absolute deviation (MAD). The dispersion 
of the various samples about the BTFR as measured by 
their MAD is shown in Fig. |4l Also shown is the me- 
dian measurement uncertainty (a). For data that scat- 



ter about the relation simply because of observational 
uncertainty, (cr) — >■ 1.48MAD as N becomes large pre- 
suming there are no systematic uncertainties that have 
not been accounted for. Fig. 2] shows that indeed (a) ~ 
1.48MAD in most cases. This leaves little, if any, room 
for intrinsic scatter in the BTFR. Indeed, the only times 
when scatter measured by the MAD significantly exceeds 
that expected from measurement errors are in the cases 
where we have identified systematic uncertainties that 
are not included in the error estimate: in the subset of 
the sample of iTrachternac h et al. (2009!) where no signa- 
ture of rotati o n has b een detected, and in the sample of 
iBegum et al.l ()2008aD where the optical and HI inclina- 
tion estimates differ greatly. Had we been unaware of 
these sources of systematic error, we might infer some 
finite intrinsic scatter. Instead, it appears that includ- 
ing less accurate data merely increases the scatter as one 
would expect. 

The BTFR of gas rich galaxies appears to have essen- 
tially zero intrinsic scatter. The same conclusion follows 
from noting that xt is close to unity (Table [2]) . A fi- 
nite intrinsic width to the relation would cause x^ > 1 
(iWeiner et al.l [2006( ) . We can place an upper limit on 
the scatter by assuming the data are Gaussian (a dubi- 
ous but conservative assumption for this purpose) and 
subtracting the measurement error in quadrature from 
1.48MAD. This gives a limit on the intrinsic scatter of 
Cint < 0.15 dex in mass for the lStarkls t al. (2009) sample 
and ajn t < 0.18 for the trimmed sample of IBegum et al.l 
(|2008al) . Essentially zero intrinsic sca tter is inferred 
by thi s method for the gas rich data of IGurovich et al. 



(2010!) or the star do minated galaxies of ISakai et ah 
gOOO). For the data of ITrachternach et al.l pOOQf ) . the 
expected error exceeds the observed scatter. This might 
indicate that the uncertainties (especially in distance) 
have been overestimated, but the sample is small so we 
should not expect the assumption of a normal distribu- 
tion to hold. Indeed, that is generally true here, which 
mitigates against there being much intrinsic scatter. 

It is common in extragalactic astronomy for correla- 
tions between physical quantities to exhibit lots of scat- 
ter. In the BTFR we seem to have a relation with little 
or no intrinsic scatter. There are some possible couple 
of caveats to this. First, it is conceivable that the un- 
certainties in the data have been overestimated. Second, 
selection effects might limit the sample to only those ob- 
jects that obey the BTFR. 

Overestimation of the uncertainties seems unlikely. 
That x^ ~ 1 is a sign not only that the model is a 
reasonable description of the data, but also that the un- 
certainties have been neither over nor underestimated. 
Even if the data were arbitrarily accurate, so that all 
of the observed scatter is intrinsic, the intrinsic width 
of the relation is still small: the observed 0.24 dex. If 
anything, this would appear to confirm that the disks of 
the selected galaxies are well behaved dynamically cold 
rotators. If it were otherwise, e.g., if there were sub- 
stantial non-circular motions in the disks of these galax- 
ies, then this would be a substantial source of intrinsic 
scatter and the observed scatter would be niuch l arger 
(jFranx fc de Zeeuwl [1991 ITrachternach et al.ll2009t ). 

It is more difficult to address the issue of potential se- 
lection effects. However, in the time since iZwaan et al.l 
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Table 2 

Statistics of Gas Rich Galaxy Rotation Curve Samples 



Sample 


N 




log A 
weighted 


log A 
unweighted 


"3 

skew 


04 
kurtosis 


(10 
weighted 


ao 
unweighted 


S 

ST 

STB 


28 
34 
47 


1.04 

0.92 
0.99 


1.71 
1.67 
1.69 


1.65 
1.64 
1.71 


-0.07 
-0.01 
-0.13 


-0.14 

0.01 

-0.10 


1.18 
1.30 
1.24 


1.36 
1.39 
1.18 



Note. — Sapmles: Si lStark et"an 120091) : T: ITrachternach et al.l 120091) : B: IBegum eFaLl 
Il2008al ). The units of A are M© km~ s* and those of ao are A s~^. The skew and kurtosis 
are the statistics of the orthogonal deviations 5 (equation (HJl . 
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Figure 4. The median error (triangles) and median absolute deviations (MAD) of the various samples. Open symbols refer to line-width 
measurements; solid symbols represent rotation curve measures. Lines connect samples that have been trimmed for quality as discussed 
in the text, showing both trimmed and untrimmed samples. Samples are labeled by first author with the MAD of gas dominated galaxies 
represented by circles and that of star dominated galaxies by squares. The MAD is a robust estimator of the dispersion of a sample; for a 
normal distribution a = 1.48MAD. The similarity of the observed scatter to the typical uncertainty is an indication that the BTFR has 
very little intrinsic scatter. 



(|1995!). the persistent surprise has been how tight the 
TuUy-Fisher relation remains when chahenged with data 
for new types of galaxies that we would naively expect 
to deviate from it. Since this exercise has now been per- 
formed by many different observers for many different 
types of galaxies, it seems unlikely that selection effects 
play a significant role in suppressing the scatter of the 
BTFR. We select gas rich galaxies here, but it is already 
well established that star dominated galaxie s obey the 
BTFR with little scatter (e.g.,|^ rhciicn 2001; iMcGaughl 
I2005bt iNoordermeer fc Verheiied 120071) . The only fur- 
ther selection made here is on data quality. When lower 
quality data are admitted, the scatter increases (Fig. [T] 
and 2]), just as one would expect. There is no clear evi- 
dence for intrinsic scatter beyond that induced by obser- 
vational uncertainty and the systematic effects inevitable 
in astronomical data. 

It is impossible to estimate the scatter that might 
result from including objects that are not already 
included. We have restricted our discussion to rotating 



galaxies, since these are the objects that define the 
TuUy-Fisher relation. It is not obvious that it should ap- 
ply to non-rotating galaxies, though the Faber- Jackson 
relation for elliptical galaxies provides an obvious 
analog. Dwarf spheroidal galaxies fall remarkably close 
to the BTFR with a simple estimate of the equiva- 
lent rotation velocity (Vf = y^a: iWolf et al.l l2010( ) 
th ough the ultrafaint dwa rfs fall systematically below 
it (jMcGaugh fc Woli[2010l ). Even objects that have no 
apparent business adhering to to the BTFR either do so, 
or come surprisingly close to doing so (e.g., the HI ring 
in Leo). Among rotating disk galaxies one would expect 
differences in the baryonic mass distribution to have 
some impact on the BTFR, just because V? oc GM/R. 
And yet disk dominated spirals, bulge dominated early 
types, diffuse low surface brightness galaxies, and dwarf 
irregulars all adhere to the same BTF R in spite of 
their palpably different mass d i stributions (iZwaan et al.1 
19951 ISpravberrv et all [iQQSt iTullv fc Verheiieri '1997: 
McGaugh fc de Blokl liggSaT iCourteau fc Rix .1999 : 
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lVerheiieDir2001l:lNoorderineer fc VerheiieDll2007[ ). Rather 
than finding greater scatter as we expand our samples 
beyond the Sc 1 galaxies traditionally utilized in distance 
scale studies, we find that rotating galaxies of all types 
adhere to the same BTFR. It therefore seems unlikely 
that there remains some large, unknown population of 
galaxies in rotational equilibrium that would induce 
substantial scatter in it. 

2.6.4. Skew and Kurtosis 

A persistent problem with interpreting astronomical 
data is whether gaussian statistics actually apply. One 
reason for this is that systematic errors often dominate 
random ones. We can check whether this might be an 
issue here by examining higher order statistics like the 
skew (lopsidedness) and kurtosis (pointiness) of the dis- 
tribution in the deviations 5 perpendicular to the BTFR. 
If the data are well behaved (i.e., distributed normally), 
the skew and kurtosis should be small. 

For all three combinations of samples in Table [21 both 
the skew and the kurtosis are small: \az\ and |d4| < 0.15. 
These samples are consistent with normal distributions 
as they should be if random errors dominate over system- 
atic effec ts. This i s espe ciall y true for the comb i nation 
of the iStarket all (|2009f ) and iTrachternach et all ()2009f ) 
samples, which have jaal = d4 = 0.01. 

For comparison, the sample of iGurovich et al.l ()2010f ) 
has az = 0.63 and a^ = —0.51 around their best fit. 
If we trim the three outlying galaxies discussed in ^2.5[ 
the skew is reduced (as = —0.15) but the kurtosis grows 
((i4 = —1.28). The negative kurtosis means that the dis- 
tribution is flat-topped: it does not have a well-defined 
central peak as it should if the data simply scatter around 
the correct relation. These data are not distributed nor- 
mally. This undermines the basic assumption of the 
least square f itting method, so the apparently straight- 
forward fit of IGurovich et al.l ()2010t ) yields a misleading 
constraint on the slope. 

3. DISCUSSION 

Understanding the physical basis of the BTFR is cen- 
tral to one of the most fundamental issues presently con- 
fronting extragalactic astronomy. The BTFR provides 
a simple test that i n principle can distinguish between 
ACDM and MOND (jMilgromll 19831) . Here we discuss the 
interpretation of the BTRF in terms of each paradigm in 
turn. 

3.1. KCDM 

The interpretation of the BTFR in ACDM depends on 
how we relate dark and baryonic mass. Different assump- 
tions can plausibly be made, leading to rather different 
interpretations. Here we start from general considera- 
tions, later proceeding to specific models. 

It is conventional in cosmology to refer to structures 
by the density contrast they represent with respect to 
the critical density of the universe. The mass enclosed 
within a radius encompassing the over-density A is 



A/a 



47rA 



-PcritRA- 



With pc 



3Hq/8ttG, this becomes 



(10) 



(11) 



By the same token, the circular velocity of a tracer par- 
ticle at Ra is V^ — GMa/Ra- Consequently, 



Ma = iA/2)-^^\GHo 



vl 



(12) 



In ACDM, the density contra st A fti 100 mar ks the virial 
extent of a dark matter halo ()Eke et al.ll2001|) . For Hq = 
72kms-iMpc^\ 



M^,r = (4.6 X 10^ Mokm^-'s-^)!^, 



(13) 



This includes all mass, dark and baryonic, that reside 
within the radius Rvir- 

Dark matter halos thus obey a scaling rela- 
tion reminiscent of the TuUy-Fisher relation (e.g., 
iSteinmetz fc Navarrolll999l ). This is not the same as the 
observed Tully-Fisher relation, as both mass and veloc- 
ity refer to a radius Rvir that is much larger than what 
is observed. In order to map these theoretical quantities 
into the observed baryonic mass Mb and circular velocity 
Vf, we need to introduce some factors: 

Mb = fdfbMv^r (14) 

and 

Vf = fvVvrr. (15) 

The total mass of baryons available within i?„i,. is 
fbM yir, where /& is the cosmic baryon fraction (/{, — 
0.17: lKomatsu et al.ll2009l ). The factor fd represents the 
fraction of these available baryons that contribute to the 
observed baryonic mass. Similarly, the observed velocity 
need not be identical to Vmr] hence the factor /y. 

3.1.1. General Constraints 

The BTFR would follow from equation (flSll with a 
slope of a; = 3 and a normalization log A = 4.74 if 
/d = /v = 1 • This is shown as the dashed line in Fig. [5j 
Had this matched the d ata, it would provide a very satis- 
factory interpretation (jSteinmetz fc Navarrol ll999). In- 
stead, it exceeds the data, leading us to consider models 
with /d < 1 and /v 7^ 1. 

Equations I13II15I to provide a generic constraint on 
ACDM galaxy formation models. These typically invoke 
feedback to suppress the cold baryon content of disks so 
that /d < 1. We define a feedback efficacy parameter 



log£ = 31og/y - log fd- 



(16) 



When fy = 1, £ is just the inverse of fd, so a galaxy with 
f = 10 has lost 90% of its baryons, and one with £ = 100 
has lost 99% of them. By "lost" we mean that they are 
not part of the disk of stars and cold gas that contribute 
to Mb in the BTFR. It is irrelevant whether these baryons 
have been physically expelled as envisioned in many feed- 
back models, or remain in the halo in some dark, unob- 
served form. £ merely tells us how effective feedback 
(or any other mechanism) has been in suppressing the 
inclusion of the available baryons into the observed disk. 
In addition to variations in/^, the mapping of the ob- 
served to virial velocity fv can also play a role. This 
ratio can easily vary with mass as more massive disks 
boost fv by engendering more adiabatic contraction in 
their host dark ma tter halos than low mass disks (e.g., 
iBullock et al.ll200ll ). My own experience with mo deling 
adiabatic contraction (jSellwood fc McGaughll2005f) leads 
me to expect this to be a modest effect, most likely in 
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Figure 5. The BTFR in ACDM. Symbols as per Fig. [T] The dashed line illustrates the ba ryons that are av ailable within the virial radii 
of dark matter halos with V-uir = Vj. Solid lines in the top left panel illustrate the model of I Mo et al.l 1)19981) for two choices of their disk 
fraction parameter rrij^. I n the top right pane l, the curved solid line is the model of Truiillo-Gonicz et al. (201(1) while the box illustrates 
the range of the model of[Tonini et al.l 1)20101). The bottom left panel shows the model of dc Rossi et al.. 1,201(1). The bottom right panel 
plots the model galaxies of Piontek & Steinmetz (2011) and Governato ct al. (2010) along with the data. Here the peak rotation velocity 
Vp is used, as this is more directly comparable to the models of iPiontek fc Steinmetz (201ll ). The inset shows that these agree quite well 
with the data for star dominated galaxies. 
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the ra nge 1 < fv < 1-3. Very recently, iReves et alj 
(|2011f) have combined dynamical and gravitational lens- 
ing measurements to estimate /y « 1.3 for galaxies with 
M^, > 6 X 10^ Mq. This represents the upper end of 
the BTFR (roughly Vf > 100 kms"!). In order to help 
with the observed variation in £", fy would need to be 
substantially less in lower mass galaxies. While these 
often have rising rotation curves (so possibly fv < 1), 
flatness of t he rotation curve i s an explicit selection cri- 
terion of the lStark et al.l ()2009l ) sample. For these galax- 
ies, /y ~ 1 seems the best guess, as their disks are not 
massive enough to substantially alter their dark matter 
halos to yield /y > 1, nor does it seem likely that the 
rotation curve will resume rising after being observed to 
flatten out so as to yield fy < 1. For the very low mass 
galaxies of lBegum et al.l (j2008aD , /y < 1 is certainly pos- 
sible, but it seems unlikely to be tiny. If fv — 0.9 at 
Vf = 30 kms~^ and 1.3 at 300 kms~^, the total varia- 
tion in f is a factor of 3. This does not, by itself, suffice 
to explain the data. 

Equation 1161 maps between the parameters of the dark 
matter halo and those of the disk galaxy it hosts. It is 
completely general, encompassing what needs to happen 
to connect theory with observational reality. In order to 
be consistent with the data presented here, models must 
meet the requirement 



log£ = 1.2 -log 



Vf 



lOOkms-i 



log 



A 
100 



(17) 



This constraint applies to models of disk galaxies over 
the range 30 < V/ < 300 kms~^. In applying this con- 
straint, care should be taken to use consistent definitions 
of total mass (Ma with A = 100) and rotation velocity. 
The term involving A is a simple shift in normalization 
that translates to other common definitions of halo mass 
(e.g., A — 200 instead of 100). The constraint holds ir- 
respective of the particular choice of definition for halo 
mass. 

Equation [T7] is a scaling between £ and rotation ve- 
locity. That the efficacy of feedback scales with veloc- 
ity is not surprising, as the shallower potential wells 
of low mass halos should be less effective at retaining 
baryons. The efficacy of feedback should also depend on 
the physics driving it. For low mass systems like those 
under consideration here, feedback is often imagined to 
be driven by star formation and subsequent supernova 
events. In this case, we might expect some correlation 
between £ and measures of star formation. 

Fig. [S] shows £ as a function of past average star for- 
mation rate and gas mass fraction. The gas fraction 
fg = Mg/Mt, while to obtain the past average star forma- 
tion rate we assume all galaxies have an age of 12 Gyr 
so that (SFR) = M*/(12 Gyr). Both of these quanti- 
ties correlate with £, albeit with substantial scatter at 
low (SFR) and high fg. The curious thing about these 
correlations is that feedback has been most effective in 
removing or suppressing baryons in galaxies that have 
experienced the least star formation and have the most 
cold gas remaining. This would seem to be the opposite 
of what one would expect from a picture in which intense 
star formation results in winds that drive out baryons. 

The effect is not subtle. Galaxies with fg « 0.9 have 
£ « 25; those with (SFR) Ki IQ-^ Mq yr"! have £ « 40. 



It appears that the vast majority of baryons associated 
with the halos of these galaxies have been ejected or oth- 
erwise suppressed. 

A more remarkable thing is that £ scales so precisely 
with Vf. In order not to impart scatter to the BTFR, 
there must be very little scatter in equation 1171 I would 
expect feedback from supernovae to be a chaotic process: 
the lack of scatter in £ seems unnatural. A halo of a 
given circular velocity could easily host galaxies of very 
different £. Apparently, they do not — all halos of the 
same Vf have basically the same £. 

The tight relation between Vf and £ provides strong 
constraints on models of galaxy formation. Indeed, it 
seems to me that it constitutes a serious fine-tuning prob- 
lem. One can readily envision models tha t obtain equa- 
tion n7\ as a mean scaling relation (e.g., iStringer et al.l 
120 llh . However, doing so with so little scatter is a much 
tougher requirement. I would ex pect fd and fv to fill out 
some distribution o f values (e.g. JKauffmann et al.lll999l : 
iBullock et al.|[20QTI ). even for halos of fixed mass. There 
is no obvious reason why these distributions should be 
either narrow or closely coupled. Yet they must be in 
order to avoid introducing substantial scatter into the 
BTFR. 

Stochastic feedback is not the only expected source of 
scatter in the BTFR. The scatter in the concentration- 
circular velocity relation for dark matter halo s is another. 
Without some tuning (Bul lock et al.l [20011 ) . this alone 
exceeds the scatter allowed by the error budget. 

Another expected source of scatter is the triaxi ality 
of dark matter halos (e.g., iKazantzidis et al.ll2010[ ). If 
the gravitational potential is non-symmetric, then the 
observed velocity depends on the viewing angle with re- 
spect to the principle axis of the potential. Otherwise 
identical galaxies will appear to have difference Vf de- 
pending on how they happen to be oriented on the sky. 
A signi ficant degree of non-axisymmetry is clearly ex- 
cluded (I Andersen et ail 1200 It Im chternach etldl 120091: 
iKuzio de Narav et al.ll2009[ ). This constraint might be 
lessened by the rounding of dark matter halos by the 
formation of the baryonic disk, but the disk needs to 
be he avy enough to have an impact. IKazantzidis et aLl 
(|2010| ) find that disks need to be at least half maximal to 
have a significant impact on their halos. The low surface 
brightness disks that make up the majority of the gas 
rich galaxy sample typically fall a factor of two below 
this level ;McGaugh (2005b). 

There are other factors though could, and proba- 
bly should, contribute to the scatt er (jEisenstein fc Loebl 
fTgM IMcGaugh fc de Blok|[T998ah . It is difficult to ac- 
commodate all the expected sources of scatter in the er- 
ror budget. This holds even if we have overestimated the 
observational uncertainties, simply because the observed 
scatter is small. 

The small intrinsic scatter in the BTFR ap- 
pears to be a serious challenge for ACDM. This is 
hardly the only p r oblem faced by the theory (e.g. , 
Kroupa et al. '20Tft_| lOliTo de Narav fc SpekkensI I2OIII: 
Bovill fc Ricotti 201l|), so the persistent failure to detect 
WIMP s is worrisome. The IXENONIOO CollaboratiOT] 
()2011f) has excluded much of the mass-cross section 
para meter space where WIMPs were expected to re- 
side ()Trotta et al.l I2008D . Given the persistence of 
dynamical puzzles (e.g., IMcGaugh fc de Blok il998al : 
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iSellwood fc Kosowskvl[200ll ). one might worry that non- 
baryonic dark matter does not exist after aU. 



Figure 6. The efficacy of feedback (equation I16II as a function of past average star formation rate (left) and gas fraction (right). Symbols 
as per Fig. [T] Feedback solutions reconciling ACDM with the BTFR require galaxies that have experienced the least star formation to have 
been most effective in driving out baryons. This holds even in galaxies where there are almost no stars (fg — > 1). 

in normahzation of the relation for galaxies residing 
in sub-halos relative to that of primary galaxies (their 
Fig. 6). This shift is already apparent at relatively large 
Vf ^ 100 kms~^. There is no evidence of this in the 
data, which includes both primary galaxies and dwarfs 
that are sometimes isolated and sometimes group mem- 
bers (and hence sub-halos). For example, KK98-251 is 
a dwarf Irregular satellite of NGC 6946, but both lie on 
the same BTFR with no indication of a shift in normal- 
ization. Current data for the ultrafaint dwarf satellites 
of the Milky Way do place them below the BTFR, but 
this shift occurs at much smaller scales (Vf S, 20 kms~^ 
Walker et al. l [20091: IWolf et all \20Wt. IMcGaueh fc"VV^ 



3.1.2. Some Specific Models 

The generic constraints discussed in the previous sec- 
tion may or may not be adequately addressed by specific 
models. There are very many relevant models in the 
literature. It is well beyond the scope of this paper to 
review all of them. Nevertheless, it is useful to examine 
s everal illustrative cases. 

iMo et ail ()1998() built a model following the same logic 
outlined above. They define a parameter rrid that is the 
fraction of the total mass in the disk. This is analogous 
to fd, such that m^ — fb when fd = ^- The resulting 
Tully-Fisher relation (their equation 4) is shown as lines 
in the first panel of Fig. [5] for nid — 0.0 5 (as they sug- 
gest) and for m^ = 0.01 (as suggested bv lGurovich et al.] 
120101) . The higher value gives a tolerable approximation 
to high mass spirals while the lower value is appropriate 
to gas dominated dwarfs. A single universal value, as 
often assumed in galaxy formation modeling, does not 
work simply because the slope of the BTFR is steeper 
than the ACDM halo mass — rotation velocity relation. 

Mor e recent mod e ls inc l ude semi-analytic trea t ments 
(e.g., iTonini et al.) [Mot iTruiillo-Gomez et al.l IMol) 



2011 l 



and f ully numerical r esults ( e.g.. iMaver fc Moord 



20041: iGovernato et al.l [20071: Ide Rossi et al.l .__ _. 
Governato et al.l 120101 : iPiontek fc Steinmetzl 1201 ID . 



These are shown in the remaining panels of Fig. [Sj 
All invoke feedback in some way, but the prescription 
for feedback is specific to each so some come closer to 
matching the co nstraint on £ than others. 

The model of iTonini et al.l (|2010f ) exhibits a substan- 
tial scatter and dependence on the distribution of stellar 
mass. This is natural: the distribution of b aryonic mass 
should matter (jMcGaugh fc de Blok|[l998al) . This model 
matches the high mass end of the observed BTFR rea- 
sonably well, but over-predicts the mass in galaxies with 
Vf ^ 100 kms" -'^. Nevertheless, this model, like that of 
IMo et ai] (|1998[ ) , provides an honest example of what we 
most naturally expect from ACDM: a substantial amount 
of scatter about a relation with a slope of approximately 
3. ^.^.^.^ 

The model of ITonini et al.l (|2010D also indicates a shift 



The model of lTruiillo-Gomez et al.l (|2010D exhibits sig- 
nificant curvature. It matches the data well at Vf ~ 
100 kms^^ but falls below it at both high and low 
Vf. There is no evidence in the data for curvature in 
the BTFR of rotating galaxies. It is of course conceiv- 
able that curvature may appear beyond the bounds of 
the current data, and may already be observed if one 
considers dwarf Spheroidals along with dwarf Irregulars 
(jMcGaugh fc Wolj|2lno[) . This occurs at a significantly 
lower velocity scale (~ 20 kms~^) than in the model of 
ITruiillo-Gomez et al.l (IMol ^ 70 km s~^). 

The model of 'de Ros si et al.l (J2010D fares better, being 
in rather good agreement with the bulk of the data for 
gas rich galaxies. Like most ACDM models, its slope 
is relatively shallow (~ 3.2; their Table 3), but this is 
only problematic at the high mass end where the model 
under-predicts the mass of star dominated galaxies. This 
might be remedied by adopting a different IMF, but the 
change would have to be substantial: a factor of 3 or so 
in stellar mass. The scatter is respectably small. 

Most recently, detailed numerical simulations have 
Toduced model galaxi es consistent with the BTFR 
Governato et al.l 120101: IPiontek fc Steinmetzl 120111) . 
These are shown in the final panel of Fig. [5l Here we 
plot the velocity measured at the peak of the baryonic 
rotation curve Vp rat her than Vf because this is c loser to 
what is reported by IPiontek fc Steinmetzl ()2011|) . Even 
though there is not a great deal of difference between 
these two measures of rotation velocity, the already good 
agreement between the models and the data shown by 
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iPiontek fc Steinmetzl ()2011[ ) improves even further when 
Vp is utihzed. 

Numerical simul ations of this sort remain c omputa- 
tionally expensive. IPiontek fc Steinmetd ()2011h present 
eight model galaxi es, which at present i s a large number 
for such models Governato et alJ ()2007| ). Of these eight, 
seven are very nicely consistent with the observed BTFR. 
As they discuss, the other case is very discrepant. Obvi- 
ously, no conclusions can be drawn from one case in eight. 
But it does raise some obvious and important questions. 
What is the intrinsic scatter of the BTFR? How often we 
should expect to see grossly discrepant cases? Would we 
recognize these objects as normal disk galaxies? 

Resolution also remains an issue in numerical model- 
ing. The model galaxies of IPiontek fc Steinmetd (|201lD 
are all in t he star dominate d regim e. Only the one model 
galaxy of iGovernato et aD (2010) resides in the regime 
of gas dominated galaxies where the comparison is most 
clean. Comparison with the Millenium simulation sug- 
gests that simulated galaxies might well follow the ob- 
served slope and normalization of the BTFR down to 
low rotation velocities (White & Guo 2011, private com- 
munication). The critical issue of the intrinsic scatter 
that should be exhibited by low mass model galaxies in 
such simulations remains open. 

Detailed numerical simulations appear to have finally 
reached the point where they can produce quantitatively 
realistic model galaxies. This is clearly a great success. 
It remains to be seen if they can provide a satisfactory 
explanation for the small intrinsic scatter in the BTFR. 
If they can, then the question becomes whether feedback 
prescriptions that accomplish this feat correspond to how 
feedback works in nature. 



served size or surf ace brightness residuals from the TuUy- 
Fisher relation (lZwa,an et al.l [l 995: S pravberrv et al.l 
[1991 ICourteau fc Rixlll999l: IMcGaughi2005ali ! 

Indeed. IMilgromI (|1983D specificallv predicted that low 
surface brightness galaxies would share the same TuUy- 
Fisher relation, with the same normalization, as high sur- 
face bright ness spirals. This prediction w as subsequently 
confirmed (jMcGaugh fc de BlokllT998bl) . The data for 
gas rich galaxies now provides the opportunity to test 
another a priori prediction. 

3.2.1. A Test with No Free Parameters 

The circular velocities of the gas rich dwarfs considered 
here are measured. Their baryonic masses are observed. 
Their location on the BTFR is fixed. Consequently, the 
gas rich dwarfs provide a test of the MOND prediction 
with zero free parameters. 

The results of this test is shown in Fig. [71 Even though 
the value of ao was set by early fits to bright spirals 
(ao = 1.2 As~^: [Begeman et al.lll99l[ ). equation IT^ fits 
the data for much lower mass galaxies well. The gas rich 
dwarfs follo w the prediction of MOND with no fitting 
whatsoever (|McGaughll201lD . 

There was no guarantee that MOND would pass the 
test posed by gas rich galaxies. A TuUy-Fisher type re- 
lation is built into the MOND force law, but its slope is 
fixed and its normalization is constrained by star domi- 
nated galaxies. The data for gas rich galaxies could have 
fallen anywhere in the BTFR plane. They might have 
exhibited the slope of 3 that is more natural to ACDM. 
Instead, they fall along the one and only line permitted 
in MOND. 



3.2. MOND 

Among alter natives to dark matter, one, MOND 
(|Milgromlll983| ) , makes strong a priori predictions about 
the BTFR, which is a consequence of the form of the 
force law in MOND. This deviates from purely Newto- 
nian at small accelerations, a < ao, where oq is the one 
new parameter introduced by the theory. The value of oq 
must be co nstrained by observation (ao w 10""'^° ms~^: 
iSanders fc McGaugh 2002) but once specified is con- 
stant. 

In the deep modified regime, a ^ ao, the effective 
acceleration a -> y/gNOo, where gN is the Newtonian 
acceleration calculated for the observed baryonic mass 
in the usual way. For circular motion around a point 
source, we can equate the centripetal acceleration to this 
effective force to obtain 



aoGAf = yl 



(18) 



Since both ao and G are constants, and all mass is bary- 
onic, one recognizes the BTFR, Mb oc Vf". 

There are several consequences of equation IT^ 
(jMilgrom|[l983l) . The TuUy-Fisher relation is absolute, 
being a direct consequence of the force law. There 
should be no intrinsic scatter. The normalization is 
set by constants of nature. Contrary to the case in 
Newtonian dynamics, the circular velocity does not de- 
pend on the radial extent of the mass distribution, 
so there is no dependence on surface brightness or 
scale length. This naturally explains the lack of ob- 



3.2.2. The Value of ao 

The data for gas rich galaxies are consistent with 
existing estimates of the acceleration constant ao 
([Sanders fc McGaughll2003) . We can reverse the ques- 
tion, and ask what the best value of ao is according to 
these data. The best fit intercept of the BTFR already 
contains this information, and the corresponding value 
of ao is given in Table [21 

There are some subtle issues to consider in calibrating 
ao . This has previously been done by fits to full rotatio n 
curves (e.g., IBegeman et al.l [19911 : ISwaters et al.|[20Tol) . 
These contain more information than simply the outer 
circular velocity, but require an additional assumption. 
In particular, a form for the function that interpolates 
between the Newtonian and MOND regimes must be as- 
sumed. For an assumed form, the best fit mass-to-light 
ratio must be found simultaneously with the best fit value 
of ao. There is not much range allowed in either of these 
parameters, but there can be covariance between them, 
so it is challenging to disentangle them uniquely. Here, 
the data are not so precise, and we make no use of the 
information in the shape of the rotation curve, only its 
outer value. However, this approach has the virtue of 
making a minimum of assumptions and using mass esti- 
mates that are completely independent of the theory. 

Strictly speaking, equation [THl applies at infinite dis- 
tance from an isolated point mass. Real galaxies are not 
point masses, and Vf is measured at a finite radius. Con- 
sequently, the intercept of the empirical BTFR is related 
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Figure 7. The BTFR in MOND. S ymbols as per Fig. [T] The 
band shows the expectation of MOND fMilgrom' 1983). The width 
of the band represents the iblcr uncertainty in ag based on d etailed 
fits to the rotation curves of a subset t Begernan et aI.iil99lD of the 
star dominated galaxies (gray squares). The gas rich galaxies test 
the MOND prediction in a new regime with zero free parameters. 



to the acceleration constant of MOND by 



A^ 



X 
aoG 



(19) 



where x is a factor of order unity that accounts for the 
fi nite size and non-spherical ge ometry of galaxies. 

iMcGaugh fc de Blokl ()1998b[) estimated x = 0.76 from 
purely geometrical considerations — thin disks rotate 
faster than the equivale nt spherica l mass distribution. 
The sample discussed bv lMcGaughl (|2005b( ) has MOND 
fits performed for oq = 1.2 As~^ and a separate fit for 
the BTFR yielding A = 50 M0km"^s*. This gives an 
empirical estimate of x = 0.80 which I adopt here. The 
gas rich galaxy data thus give 



flo 



= 1.3±0.3 As-2. 



(20) 



The uncertainty here includes the formal uncertainty in 
A (equation |9]) plus an allowance for 20% uncertainty in 
X- It is conceivable that there could be a small system- 
atic offset in the mean value of x between gas and star 
dominated galaxies. 

Detailed rotation curve fits to some of the same galax- 
ies suggest ao w 0.7 As"^ (jSwaters et al.ll2010f ). This 
becomes ao ~ 1.0 As"'^ if we exercise the same qual- 
ity control on the data as applied here. This gives some 
idea of the systematic uncertainty and the impact of the 
best-fit stellar mass-to-light ratio when it is treated as a 



fitting parameter along with oq. However, data quality 
remains the paramount issue. 

The BTFR is clearly a success of MOND, despite 
its many other drawbacks. Further tests are possi- 
ble, as adherence to the BTFR does not guarantee 
that it is possible to obtain an acceptable fit to the 
entire rotation curve of each galaxy. While MOND 
is kn own to do quite w ell in some gas rich galaxies 
('e.g. JGentile et al.ll2010D. it may have problems in oth- 
ers (e.g.. lSanchez-Salcedo fc Hidalgo-Gamezil201lL where 
again the issue appears to be with the inclination). It is 
particularl y challenging to fit the full rotation cur ve of 
NGC 3198 (IBottema et al.l[2002l: [Gentile et al.l[201ll) . the 
most massive galaxy to meet the gas rich criterion im- 
posed here. Perhaps this suffices to falsify MOND, but 
we should not be so eager to disbelieve MOND as a the- 
ory that we ignore the remarkable empiri cal vi rtues of 
the simple formula suggested by Milgromi (|1983f ). 

There are other suggestions for modifying grav- 
itational theory besides MON D . Co nformal Weyl 
gravity (iMannheim fc KazanasI |1989() and MoG 
(jMoffatl l2006f ) offer two examples. If these theories 
can fit rotation cur ves ([Brownstein & Moffat 200y; 
IMannheim fc O 'Brien 2010), then they wiU also fit the 
BTFR ()Moffat fc Toth.201li However, I am not aware 
of a clear a priori prediction of this phenomenon by any 
theory other than MOND. 

4. CONCLUSIONS 

The baryonic mass-velocity relation provides an im- 
portant test of ACDM and alternatives hke MOND. I 
have assembled recent data for gas rich rotating galax- 
ies in order to test these theories. With Mg > Af,, the 
location of these objects in the baryonic mass-rotation 
velocity plane is effectively measured directly from obser- 
vations without the usual systematic uncertainty in the 
stellar mass-to-light ratio. 

A long standing prediction (|Milgroinlll983[ ) of MOND 
is that rotating galaxies will fall on a single mass-velocity 
relation with slope 4: Mb (x Vf. This prediction is real- 
ized in multiple independent data sets. Gas rich galaxies 
fall where predicted by MOND with no free parameters. 
There are not many predictions in extragalactic astron- 
omy that fare so well over a quarter century after their 
publication. 

In ACDM, the naive expectation of a of Tully-Fisher 
type relation with a slope around 3 fails to predict the 
location of gas rich galaxies in the plane of the BTFR. 
Only the m ost recent detailed numerical models (e.g., 
lGoverna"to" et al. 2010; Piontck & Stcinmetz 2011) suc- 
ceed in reproducing the observed phenomenology. This 
is certainly progress. However, a physical understanding 
for why galaxy formation in the context of ACDM should 
pick out the particular phenomenology predicted a priori 
by MOND remains wanting. 



The work of S.S.M. is supported in part by NSF grant 
AST 0908370. 
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